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Abstract 

The nuclear spin-rotation interaction in the hyperfine structure of the hydro- 
gen molecular ion is investigated. The interaction constants are determined 
and are found to differ in sign and magnitude compared to another theory, 
but they are in agreement with some values derived from experiment. 
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Information on molecular structure nuclear forces 0,^, fundamental symmetries 0,0, 
and even interstellar molecules P] can be gleaned from experimental and theoretical studies 
of molecular hyperfine structure. The precision of measurements of diatomic hyperfine 
transition frequencies is continually improving because of developments such as ion trap/rf 
spectroscopy [j^, ion beam/laser beam and laser-radiofrequency double-resonance 
methods. Moreover, a recent proposal for the trapping and cooling of paramagnetic 
neutral molecules should offer — once spectroscopy has been carried out — the ultimate in 
precision: natural linewidth resolution of hyperfine frequencies. 

For molecules, which have no net electronic spin, the interaction energy of each nu- 
clear magnetic moment with fields generated by the motion of the other charged particles 
in the molecule is important. It has been studied in detail both experimentally and theoret- 
ically for a number of such molecules (cf. Ramsey |]12[ and Townes and Schawlow The 
primary effect is the nuclear spin-rotation interaction originating with the magnetic fields 
generated by the rotating nuclei and orbiting electrons. The interaction constant can be re- 
lated to the magnetic shielding constant [ll3|,p!4| , which describes the effective magnetic field 



at a nucleus in an external magnetic field and is the basis of chemical shifts in NMR spec- 
It is also a sensitive test of electronic wave function calculations 



14,16,17 



troscopy [0,0 

Recently, nuclear spin-rotation interaction has been interpreted as a (Berry or geometric 
phase-like) manifestation of a nonabelian gauge potential in molecular physics [0. 

Nuclear spin-rotation interaction is also present in molecules, which have a net elec- 
tronic spin, but it causes a much smaller energy in comparison to effects arising from the 
interactions of the electron spin and thus it is usually not included in the phenomenological 
spin Hamiltonians used for fitting measured hfs transition frequencies. Some exceptions are 
and alkaline earth monofluorides [ 2C ,Q , where empirical values for nuclear 



H2 + 

spin 



No+ 19 



rotation interaction constants have been obtained. 
In this Letter, the nuclear spin-rotation interaction for the hydrogen molecular ion H2''" 
is studied. It is found that for the ground vibrational state the interaction constant is 
—41 kHz, nearly | the value for H2, —113 kHz [0, but differing from a previous theory 
in both magnitude and sign. For higher vibrational states there is agreement with some 
empirical data. 

Theories giving nuclear spin-rotation interaction constants for diatomics with other than 
states exist, but they are complex. A simplification for H2^ is obtained by utilizing the 
theory for H2. The isotropic magnetic shielding constant a, which describes the net magnetic 
field seen at a nucleus in an external uniform applied magnetic field, is related to the nuclear 
spin-rotation interaction constant. The shielding constant for nucleus "a" can be written 



a{R) = ai^{R) + aM{R), 
where, in atomic units, R is the internuclear distance, and 



(1) 



a^R) = ia'(0|l/r,|0) 



(2) 



is the "Lamb" [^,0 or diamagnetic part, where joins the nucleus to the electron and 
|0) is the ^E+ electronic wave function, and crhf(-R) is the high-frequency or paramagnetic 
part. The designation high-frequency arises because the expression for a^f involves 
highly-excited electronic states [O], 
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MR) 



a 
~6 



E 



(0|L,K)-(i|r-^LjO)^^^._ 



Eo{R)-E,{R) 



(3) 



where the symbol J2i^o represents an infinite summation-integration over the intermediate 
electronic states with wave functions "adj." indicates the Hermitian adjoint of the 
preceding term, Eq{R) and Ei{R) are the respective electronic energies and is the orbital 
angular momentum of the electron about the nucleus "a" , for which a is being evaluated. 

The separation (|l|) into a diamagnetic and paramagnetic part is not unique [|13],|14| . It is 
related to the choice of origin for the electron orbital angular momentum and other vectors, 
which in turn is related to the choice of a gauge constant associated with the magnetic vector 
potential in the Coulomb gauge (cf. [0,^,^). Expressions (0) and (0) above adopt the 
convention of the vector origins at the nuclei. 

Let I and K be, respectively, the total nuclear spin and rotational angular momenta 
and define the nuclear spin-rotation interaction constant /, a frequency, through the energy 
hfl ■ K in the hfs Hamiltonian P^p5[ , where h is Planck's constant. Introducing, 

(4) 



fiR) = MR) + UR), 



the major contributions [|T2|,^ to the energy are from the interaction of each nuclear mag- 
netic moment with the magnetic field generated by the other rotating nucleus. 



hfi{R) 



i?3 ■ 



and with the magnetic field generated by the orbiting electron |13], 



hh{R) 



a- 



o"hf(-R), 



(5) 



(6) 



where R is the internuclear distance, and the dimensionless quantities a, Qp ~ 5.586, and 
(Thf are, respectively, the fs constant, the proton (yf-factor, and the high-frequency component 
of the magnetic shielding constant as defined in Eq. Atomic units are used throughout, 
except for /, which is expressed in kHz. 

Direct computation of o"hf(-R) is rather involved and requires evaluation of a term in 
second-order perturbation theory. It is much easier to use (|T]) because values of a are 
available and values of ctl are easily calculable. Moreover, because a is gauge-independent 
the irreducible components of the (symmetric, second rank) shielding tensor calculated with 
a different gauge origin by Hegstrom ||2^ could be utilized here to calculate a p8|. Values 
of crL(-R) were calculated |^ at various internuclear distances R and are given in Table |. 



Then, (Xhf was determined using (|1|) and and values of / were determined using (^. 
Values of a and / are given in Table |. 

Jefferts measured hfs transition frequencies of the vibrational v = 4-8 states of the 
ground electronic state of H2'''. For the rotational quantum number K = 1 for each v, he fit 
the transition frequencies to the Hamiltonian 



Hms = bI-S + chS, + c/S ■ K + /I ■ K 



(7) 



and obtained values for the coupling constants 6, c, and /, where I and K have been 
defined, and S is the electronic spin angular momentum vector, and his results for / are 
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given in the Table |T|, col. 2. The measured hfs transition frequencies were refit to the 



Hamiltonian (J^) independently by Kalaghan Menasian and Varshalovich and San 



nikov who all obtained mutually consistent results for / that had the opposite sign 



and somewhat different magnitude than those of Jefferts. The values of are given in 
Table |II|, col. 3. 



McEachran, Cohen, and Veenstra |Q obtained a theoretical expression for a nuclear 



spin-rotation interaction constant by simply multiplying the first order electronic spin- 
rotation constant di ||3^ by the ratio of the electron mass and the nuclear reduced mass. 
Their formula gives a positive interaction constant. Interestingly, it is independent of the 
proton (yf-factor Qp, and thus does not contain the proton magnetic moment. Nevertheless, 
close agreement was obtained with the empirical values of / obtained by Jefferts. The values 
for V = 4-8 of McEachran et al. and the f = value calculated using their formula are 
given in column 5 of Table ||. 

Lundeen, Fu, and Hessels extracted hfs constants for from an analysis of their 
measured transition energies for highly-excited Rydberg states of the hydrogen molecule. 
Their value for / for the v = K = 1 state is given in column 4 of Table 0. It is consistent 
with zero. 

In order to compare the present results with the above, the values of / in kHz — a partial 
listing is given in Table |T| — were averaged over the vibrational-rotational wave functions 
for various v with K = 1 calculated with the Born-Oppenheimer potential using standard 
methods. The results are given in the last column of Table ||. Good agreement, including 
the sign, is obtained with the empirical results of [^0| -|3^ . 



The expression (^) can be obtained using more formal arguments. The first term in (^, 
/i, follows by consideration of the interaction of a nuclear moment with the magnetic field 
generated by the rotation of the other nucleus [^]. The second term in (^, /2, would be 
expected to occur in second order perturbation theory through the electron orbital-rotation 
interaction 

Hi = • K (8) 

where L is the electronic angular momentum about the center of nuclear mass and Mp is 
the proton mass, and the electron orbital-nuclear spin interaction 

H[ = HI + H[,, (9) 

where 

H'la = 2gpii^HNla ■ r'^^a (10) 

and 

= 2gpfiBfiNlb ■ '^fe'^Lfe, (11) 
where and lb are the nuclear spin angular momenta, and where the terms H4 and H[ were 



derived by Dalgarno, Patterson, and Somerville |^ from a non-relativistic reduction of the 
Dirac eq. for H2^. Writing 
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L,) ■ K = i^4a + 



Ahi 



(12) 



where is the electronic angular momentum about nucleus a, and similarly for L^, the 
corresponding energy in perturbation theory is, schematically, T = + T^, where 



Eo{R)-E,{R) 



+ adj. 



(13) 



and similarly for T;,. It follows that 



T 



i?2 



E 



(0|Ljz) ■ (^Ir-^LJO) 



+ adj. 



la K, 



(14) 



la + I„ and T = hhl ■ K, 



the desired result Eq. (H) is 
although a more 



and similarly for T^,. Since I 

obtained by substituting Eq. (||) for dhf into T^, Eq. (|l4|), and into T^, 
careful derivation is desirable. 

Elaborate treatments of the Hamiltonians for diatomic molecules in other than 
states including all angular momenta have yielded nuclear spin-rotation interactions for 
various electronic states 



37-41 



In particular, terms similar to /i and /2 were given by 
Mizushima |Q for a state, although they were not related to Ramsey's theory. 

The present approach may be generalizable to other molecules. Hegstrom has 
shown for that there is a relation between the second order electron spin-rotation in- 
teraction constant d2 and a second order part of the shielding constant a. Eq. (H) would 
thus imply that there is a relation connecting c/2 and /2. Physically, such a relation might 
be anticipated because the nuclear spin and electron spin are both separately coupled to 
rotation through the excitation of electronic angular momenta. Such a relation would make 
it possible to estimate /, if d were measured, using theoretical values of /i and di. Moreover, 
through the above arguments, the recent findings that electron spin-rotation in a param- 
agnetic molecule and nuclear spin-rotation in a diamagnetic molecule are each, separately, 
describable using a nonabelian gauge potential method in molecular physics |^ , p!8[| can 
probably be related. 

The present paper demonstrates that there is a solid theoretical basis for nuclear spin- 
rotation in a molecule. A synthesis of the theories for molecules and ^E molecules has 
yielded a simple expression for the interaction constant in H2'''. Using calculated properties 
of H2''' numerical values obtained for the interaction constant were found to be in good 
agreement with empirical results. The sign of the effect is found to be the same as for 
II2, making H2^ another rare example of a diatomic with negative nuclear spin-rotation 
interaction. 

The author is grateful to Prof. K. Burnett and Prof. A. Dalgarno for comments on the 
manuscript. This work was supported in part by the National Science Foundation through 
a grant for the Institute for Theoretical Atomic and Molecular Physics at the Smithsonian 
Astrophysical Observatory and Harvard University. 
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TABLES 



TABLE L Electronic expectation values of the shielding constants (Th{R) and cr{R), dimension- 
less, and values of f{R), in kHz. 



I ) I 

R/ao 


cri^{R) X 10^ 


cr{R) X 10^ 


fi.R) 


1.0 


21.5183 


17.6401 


-453.7 




1Q 34Sfi 


1 5 3584 


—91 3 8 


1.50 


17.6258 


13.6704 


-114.5 


1.75 


16.2484 


12.3995 


-67.21 


2.0 


15.1295 


11.4296 


-42.31 


2.25 


14.2116 


10.6825 


-28.16 


2.50 


13.4530 


10.1046 


-19.63 


2.75 


12.8232 


9.6575 


-14.22 


3.0 


12.2991 


9.3133 


-10.65 


4.0 


10.9538 


8.6124 


-4.285 


5.0 


10.3478 


8.4913 


-2.215 



TABLE IL Comparison of experimental and theoretical values of the nuclear spin-rotation 
coupling constant for for various vibrational states v with rotational quantum number K = 1, 
in kHz. The column labeled "Refit" gives results obtained by refitting the raw experimental 
transition frequencies from Q. Numbers in parenthesis are quoted experimental uncertainties. 



Empirical 



Theory 



Expt. II 



Refit^ 



Expt. 13 
-3(15) 



M 

46 



Present 
-41.7 



4 38(1.5) -34(1.5) 38.8 -36.2 

5 36(1.5) -33(1.5) 36.2 -34.6 

6 34(1.5) -31(1.5) 33.6 -32.9 

7 32(1.5) -29(1.5) 31.1 -31.0 

8 30(1.5) -27(1.5) 28.6 -29.1 



^Values are from Ref. |30| and are identical with those given in Refs. |31| and |32|, except |31| gives 
-32(1.5) for v = 5. 

^Values for v = 1, 2, and 3 are, respectively, —40.8, —39.4, and —38.0. 
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